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Abstract 

K. Godel [Gl] discovered his celebrated solution to Einstein equations in 1949. 
Additional contributions were made by Kundt [K] and Hawking-Ellis ([H-E], 5.7). 
On the other hand, a general Lorentz invariant operator, associated to the so-called 
"relativistic diffusion" , and making sense in any Lorentz manifold, was introduced by 
Franchi-Le Jan in [F-LJ]. Here is proposed a first study of the relativistic diffusion 
in the framework of Godel's universe, which contains matter. 

1 Introduction 

K. Godel [Gl] published his celebrated exact solution to Einstein equations in 1949. 
The most striking feature of this cosmological model was to be non-causal (though locally 
and geodesically causal), containing closed timelike curves. For this reason, it is generally 
considered as rather unphysical. Possessing a series of interesting properties, it aroused 
however a great interest among physicists. For example, it contains rotating matter, but no 
singularity. Moreover, the explicit exact solutions to Einstein equations are not so many. 

W. Kundt [K] studied its geodesies, and S. Hawking and G. Ellis ([H-E], section 5.7) 
stress on coordinates (defined by Godel himself) showing up its rotational symmetry (about 
any point), to draw a nice picture of its dynamics. D. Malament ([M1],[M2]) calculated the 
minimal energy of a closed timelike curve. K. Godel [G2] discussed other rotating universes, 
which are spatially homogeneous, finite, and expanding, and he showed in particular that 
there exist a lot of strongly causal such cosmological models. 

The purpose of the present work is to study, in the framework of Godel's universe, the 
behaviour, and mainly the asymptotic behaviour, of the so-called relativistic diffusion. 

The relativistic diffusion was introduced by J. Franchi and Y. Le Jan in [F-LJ], in the 
framework of general relativity, on an arbitrary Lorentz manifold, as the only diffusion 
which is invariant under Lorentz isometries. In this sense, it is the Lorentzian analogue of 
the Brownian motion on a Riemannian manifold. It lives in fact on the pseudo-unit tangent 
bundle of the considered Lorentz manifold, and is roughly the integral of Brownian motion 
of the unit pseudo-sphere of the tangent space. It can also be seen as a random perturbation 
of the timelike geodesic flow. 
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This article begins with a detailed study of timelike and lightlike geodesies, in a different 
and more complete way as Kundt [K] did. As a conclusion of the study of lightlike geodesies, 
a definition (Definition [2]) of light ray (or boundary point, as an equivalence class of lightlike 
geodesies, without use of causality) and of convergence to a light ray is given, which appears 
to be rather natural in this non-causal universe (it can be reinforced to a certain extent : see 
Remark [TO]) . Thus the set of light rays has a natural structure of 3-dimensional boundary, 
on which the isometry group of Godel's universe does operate. 

Then the relativistic diffusion of Godel's universe is introduced. In order to study such 
a 7-dimensional diffusion, some sub-diffusions are considered, of dimensions 1, 2, and 4. A 
leading concern is here to bring out all asymptotic variables of the relativistic diffusion, or 
in other words, the tail cx-field of its natural filtration (which is in turn closely related to 
the Poisson boundary of the relativistic diffusion). The clue in this direction is the idea, 
suggested and partially worked out in [F-LJ] in the case of Schwarzschild solution, and very 
recently established in [B] in the flat case of Minkowski space, that convergence to a light 
ray should eventually occur. The following statement, progressively established below, 
shows that this general guess stands out as reinforced, also in the case of a non-empty (and 
even non causal) universe. 

The main results of the present article are summarised in the following. 

THEOREM (i) The relativistic diffusion is irreducible (on its 7-dimensional phase-space). 

(ii) Almost surely, the relativistic diffusion path possesses a 3-dimensional asymptotic 
random variable, and converges to a light ray (in the sense of Definition [2] and Remark [TOT) . 

(Hi) The support of possible light rays the relativistic diffusion can converge to, is the 
whole 3-dimensional boundary space of light rays. 

As a consequence of this theorem, and on the basis of some secondary results and 
considerations, the following conjecture appears to hold likely : by the showing up of the 
3-dimensional asymptotic random variable evoked in (ii) above, the whole tail cr-field of 
the relativistic diffusion of Godel's universe, and then its whole Poisson boundary, has been 
brought out. 

The only relativistic case in which the analogous statement has been proved, up to now, 
is Minkowski space, in [B], by two different methods: Doob's /i-process conditioning and 
then couplings, making use of an explicit expression of the laws of already found asymptotic 
variables ; or alternatively : study of the random walk associated with a lifted relativistic 
diffusion on some (Poincare) fixed locally compact group. The use of both methods does 
not seem to be easy in the present curved case (likely as in any other curved case), since 
neither are explicit the laws of the asymptotic variables, nor appears any Poincare-like 
group. 
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2 Godel's pseudo-metric 

Definition 1 The Godel's universe is the manifold R 4 7 endowed with coordinates £ := 
(t, x, y, z), and with the pseudo-metric (having signature (+, — , — , — ) ) defined by : 

ds > ;= dt 2 -dx 2 + \ e 2 ^ x dy 2 + 2e^ x dtdy-dz 2 , 

for some strictly positive constant uj . 

The inverse matrix of this pseudo-metric ((py)) is as follows : 



((g ij D 



( -1 2e~ V2ujx \ 
0-10 

2e -V2u>x o -2e" 2 ^^ 

V -I) 



The unit pseudo-norm relation, defining proper time s , is : 



(0) 



1 + i 2 + x 2 + z 2 



e V2. Xs y s + 2 i s 



The isometry group of Godel's universe is the five-dimensional Lie group generated by : 

1) the translations (t, x, y, z) \— > (t + t , x,y + y , z + z ) of the linear (t, y, z) 3-subspace ; 

2) the hyperbolic dilatations (t, x, y, z) i— > (t,x + x ,y e'^ 2 ^^ , z) ; 

3) the rotational symmetries (u, r, 0, z) i— > (u, r, <p + O ? z), in the new coordinates system 
(u, r, 0, z) E R x R + x (R/^Z) x R defined by |i - u| < tt/cu and : 



e V2u>x = c h(2r)+sh(2r)cos(w0); e j2 " x uy = sh(2r) sin(cu0); tg [fi^(0+i-u)] = e~ 2r tg p 



2 J ' 



we have indeed ds 2 = [du + 2 sh 2 r d0] 2 - 2w 2 dr 2 - | sh 2 (2r)d(j) 2 - dz 2 . 

Godel ([Gl], section 4) proved that these three types of isometries generate indeed the 
full isometry group. As the action of this group is clearly transitive on R 4 , Godel's universe 
is an homogeneous space-time. 

Letting uj go to 0, we recover Minkowski space as limit of Godel's universe. 



2.1 Timelike geodesies 

Geodesies are associated with the Lagrangian L(£, £), given by: 



2L(i s ,Q=t 2 s -x 2 s + e 



1 a 2V2u)x 3 -2 



The equation of geodesies 



d r dL(t,Csh dL(t,i 



ds 



(1) i s + e 



a ; (2) e 



2V2w; 



y s + 2e 



t s y s - z s 



reads here : 



b ; (3) z s 
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(4) x s + (lo 
for constant a, b, c. 

Equations (1) and (2) jointly are equivalent to: 

(1') i s = be-^ Xs -a ; (2') y s = 2ae'^ Is - be~ 2V * ux ' ; 

and then using Equations (1'), (2'), (3), we see that Equation (0) is equivalent to: 

1 2 



1 + 



he -V2^x 3 _ a 



+ x 2 + c 2 = \b 2 e~ 2 ^' 



or equivalently : 



(0') x* 9 + \\2a-be-^ ux * 



J s 1 2 
,2^1 , „2 



2 2 1 

a — c — 1 



Note that necessarily a > 1 + c , and a b > . 
Then, owing to Equation (2), Equation (4) is equivalent to: 



(4') 



V2 
uj b 



x + y = Y , for some constant Y. 



Note that Equations (1) and (2) imply: 

1") i s + V2ue V2uJXs x s y s + 2V2uji s x s 



0; (2") y s -2^2ue-^ x H s x s = 0, 



so that, using the derivative of unit pseudo-norm Relation (0), we see that Equation (4) is 
implied by Equations (0), (1), (2), (3), unless x s = . 

Setting R := ajl - (1 + c 2 )/a 2 and k := R/(V2a) G [0, ^-], we must have by (0') : 

x s = R cos(oj (f s ) , 5e _v/2w:Cs = 2 a — \[2R sin(uj ip s ) , 
for some angular component ip s , whence : 

(fis = be~^ ux ' = 2a- V2Rsm(iu<p s ), 

and then : 



2au (s - s ) = J — 



ui dip 



k sin (a; ip) y/l - k 2 



Arctg 



tg(ujip s /2) - k 



VT^k 2 



Therefore 
and 

(5) 



tg (u = vl — k 2 tg a vl — k 2 uj (s — s ) +k, 



-V2 



2k(Vl-k 2 tg 


a Vl — & 2 w (s — so) 




i + (vi - & 2 tg 


a \Jl — k 2 ui (s — s ) 


+ A;) 2 



or equivalently : 

(5') 



.)■., = —J— log 



V2 



2a 



(l-k 2 ) (l+tg 2 


a y/l — k 2 ui (s — s ) 


) 


1 + (VI - k 2 tg 


a y/1 — k 2 cu (s — so) 


+ k) 



Moreover by (1') we have : 

i s = a — \Pl R sin(u v? s ) 

whence 



- a, 



2 

(6) t s = T - a(s - s ) + - Arctg(Vl - k 2 tg [a Vl - k 2 lo (s - s )l + fc). 
In this formula (6), the successive determinations of Arctg, at the successive values s E 



so + 



Z , are understood to be chosen conveniently, in order that the absolute time 



coordinate (t s ) be continuous, as it must be. Observe that (t s ) is strictly monotonic if and 
only if k < \ , or equivalently, if and only if (1 + c 2 ) < a 2 < 2(1 + c 2 ). 



Finally, by (4') and (5') we have : 

2ak 



(7) 



y s = Y + 



Aak/(bu) 



1 + (yi - k 2 tg [a VI - k 2 uo (s - s )] + k) 



2 ■ 



a 2 k 



which is consistent with (2') : y s = 4 —j— [1 — k sin(uj ip s )] sin(cj <p s ). 

Observe from Equations (5) and (7) that every timelike geodesic has a bounded, periodic 
projection in the (x, y)-plane, and moreover that it obeys the following relation : 



(8) 



-V2l 



2a 



1 



+ 



to b 
2~a 



(Vs ~ Y) 



= k . 



Remark 1 The case k = is particular. It implies (using Equations (5), (0'), (1'), (2')) : 
i 2 = 1 + z 2 and x — y — , and then : 

(x s , y s ) constant and t s = t + a s , z s = z + c s , with a 2 = 1 + c 2 . 

Reciprocally, if x = y = , then (by Equations (2'), (0')) the corresponding geodesic 
must satisfy also k — , and then be included in the phase subspace defined by : 

S ■= {x = y = 0} = [i 2 = 1 + z 2 ; x = y = o}. 

Therefore, the case k = corresponds to the geodesically stable phase subspace £ . 

Remark 2 Every timelike geodesic is defined for all proper times s , unbounded and 
causal. Moreover, it never accumulates near its past. Indeed, if for proper times s < s' 
we had x s i = x s , then by Equation (5) we should have : 



tg aVT — k 2 uj{s' — so) = tg ay/l — k 2 uj(s — s ) , then s' = s + 
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with n E N*, 



whence by Equation (6) : t s > —t s = sign(a) y^j 1 — ^Jr-k? ) ' anc ^ then 

\t s ,-t s \ = — (2 - (1 - k 2 )- 1 / 2 ) > -. 
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The following statement, which will be used later to ensure the irreducibility of the rela- 
tivists diffusion, shows up the non-causal structure of Godel's universe, despite the preced- 
ing remark[2]: the causal past of any point of Godel's universe is the whole Godel's universe. 
In particular, the causal boundary, in the sense of Penrose (or Geroch-Kronheimer- Penrose, 
see ([H-E], section 6.8)), reduces to a single point. 

Proposition 1 The Godel's universe (iR 4 , ((<7ij))J is geodesically transitive: any two 
points of it can be linked by a piece-wise geodesic timelike continuous path. 

Proof Observe from Remark [1] that (taking k = c = 0) there are timelike geodesies 
moving at will the coordinate t , without changing any other coordinate, and that (taking 
= 0, c / 0) there are timelike geodesies moving at will the coordinate z , without 
changing the coordinates (x,y). 

Hence it will be sufficient to move piece- wise geodesically the coordinates (x, y) , forget- 
ting henceforth the coordinates (t,z). 

Observe that any given geodesic can move the coordinate x only by a uniformly 
bounded value, since by Equation (5) we have : 



e 



-V2ui (x s -x s i) 



[1 + tg 2 {uj <p s )] (l + [VT^W tg (u yy) + k} 2 ) 

(l + [VT^¥ t g (u Vs ) + W) [l + t g v w)] 



so that several geodesic arcs are needed. Fix coordinates (x,y), which we want to move 
to other fixed coordinates say (x', y'). By a finite number of geodesic moves, according to 
the equation displayed just above, we can move (x,y) to (x',y"), for some y". 

Then, the quotient in the equation displayed just above equals 1 as soon as 

ip s + <f s > = uj' 1 Arctg \(2kVl - k 2 )/(l + k 2 - Vl - k 2 

so that we can choose a geodesic arc, independently of the value of the parameter b , such 
that this holds and such that [y/1 — k 2 tg (u </v) + k] 2 [Vl — k 2 tg (u (p s ) + k} 2 . Finally, 
on such geodesic, by Equation (7) we have : 



4ak 

y s - y s > = -r — x 



1 



boo 11+ [VT^k 2 tg(uip 8 ,) + k} 2 1 + [y/T^k 2 tg (uj (p s ) + k} 2 
proving that, choosing conveniently the parameter b, we can move (x',y") to (x',y' 



6 



2.2 Lightlike geodesies 

Equations (1), (2), (3), (4) remain the same, while the pseudo-norm equation (0) is re- 
placed by : 

(0") t 2 s + x 2 s + z 2 s = \ Y~^ x ' y s + 2 i s ] \ 
As previously, knowing Equations (1) and (2), Equation (4) is equivalent to 

y/2 

(4') — - x + y — Y , for some constant Y, 
ujo 

and is implied by Equations (0"), (1), (2), (3). Thus lightlike geodesies are the solutions to 
the system : 

(1') i a = be~ V2 ' ux ' -a; (2') y s = 2 a e -^ ux ' - be~ 2V * ux ' ; (3) z s = c ; 

and 



(0'") x 2 s + \\2a-be-^ x 



2 2 

a — c . 



Hence we must have again a 2 > c 2 , and ab > (if not, the trajectory must be constant), 
and, setting now 



I(l- C 2/ a 2) e[0 ,-y, 

we get the same equations (5), (6), (7), (8) as above, merely with k replacing k . 

As the parameter s cannot any longer have here a meaning like proper time, but stands 
only for an affine parameter, determined up to a change s i— > u s + v , the constant (a, b, c) 
is now irrelevant. Note that on the contrary, the constant Y of Equation (4') is consistent. 
The only meaningful a priori parameter for a lightlike geodesic is the "impact" parameter : 

(9) B = (I, q, Y) := (-, -,Y)eB := [-1, 1] x R* + x R . 

Eliminating s , we see indeed that a lightlike geodesic having impact parameter B solves : 

£dt = (ge- V2ux -l)dz; {2- ge- V2ux )e- V2ux dt = {Qe~^ x - 1) dy ; 

( ee -^ ux -l)dx = ±^l-£ 2 - \{2- Q e-^<»*) 2 dt. 

Let us sum up the description of lightlike geodesies in the following statement. 

Proposition 2 Any lightlike geodesic (x T ,y T , z T ,t T ) having impact parameter 

B = (£, g,Y) e B satisfies, for an additional parameter (Z , T ) E IR 2 and for any real r : 



e = £ x ( ! _ 2v^(^rT^tgr + ^T^) x _ 
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2 
Q 



V2(l-l 2 ) / 



Y + 

LOQ 



2 + (y/m 1 tgr+^r 

^0 



2 ' 



r 2 



t T = T +-Arctg V(l + ^)/2 tgr+ V(l-^ 2 )/2 ; 

Wv /(l+£ 2 )/2 ^ 



C 



B 



^„-^*r_ll 2 ^ — 12 1-/2 

2 



+ 



U P ( v\ 



Remark 3 The last equation in Proposition [2] shows that to any given lightlike geodesic 
is associated a cylinder Cb, parallel to the (t, z)-coordinate plane. Reciprocally, by Propo- 
sition [2] again, any lightlike geodesic which is drawn on the cylinder Cb has a prescribed 
projection on the (x, ^-coordinate plane (up to changing affine parameter r). The equa- 
tions displayed in Proposition [2] define a lightlike geodesic associated to any given B G B. 

Considering then any continuous angular parameter (p — (p T (determined modulo 2h/uj) 
such that : 

tg(w^/2) = y/(l + P)/2 tgT + y/(l-P)/2, 
then by Proposition [2] we have : 

ge -y/2u,x T = 2 _ v / 2 (i_ £2) sin^^) anc l uj q y T = uj qY — ^2(1 — £ 2 ) cos(tu^ T ), 

together with : 

t T = T . T + y? T , z T + ^t T = Z + ^T + ^ T . 

0^(1 + ^ 2 )/2 

Since the function r i— ► tp T — 2 r is 7r-periodic, the functions 

r^ut T -2{l-[2{l + f)}- l l 2 )r and r ^ z T - t Ht = Z - ^^f^ 

1 L 1 JJ ^ ^2(1 + *) -1 v^ 3 )- 1 

are 7r-periodic too. This implies that t T wanders out to infinity, nearly linearly, and that 
the projection on the (t, z)-coordinate plane of any lightlike geodesic has an asymptotic 
direction : 

hm — = , = . 

T ^ ±oc U yj2(l + I 2 ) - 1 

This prescribes geometrically the sign of parameter i , which was not determined by the 
cylinder Cb alone, which however determines Y). Note that the impact parameter 

B = (£, g,Y) £ B has thus indeed a clear geometrical meaning. 

Note finally that the additional parameter (Z ,T ) depends on a translation on the 
parameter (r,ip T ), and then, contrary to B — (£, g,Y), is geometrically irrelevant. 
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Recall that in a strongly causal space-time, it seems natural to use the causal boundary, 
in the sense of Penrose, to classify lightlike geodesies by gathering in an equivalence class, 
called a light ray, all geodesies which converge to a given causal boundary point (having 
asymptotically the same past, see ([H-E], section 6.8)). On the contrary, in the present 
setting (recall Proposition [T]), such classification is totally inoperative. It seems that no 
alternative classification has been proposed so far, which be relevant in a non-causal setting. 

However, owing to the above remark[31 we are led to adopt here the following alternative 
classification of lightlike geodesies into light rays, and then also, to see the 3-dimensional 
space of light rays as an alternative notion of (non-causal) boundary, as follows. 

Definition 2 Let us call light ray , or boundary point , of Godel's universe, any equivalence 
class of lightlike geodesies, identifying those which have the same impact parameter B = 
(£, g,Y) G B . Thus B = [—1,1] x W* + x R is the boundary of Godel's universe. 

Let us say that a path s i— > £ s = (t s ,x s ,y s , z s ) of class C 1 in Godel's universe 

converges to the light ray B = (£, g, Y) if, setting : 

a 8 := i a + e^ u *-y a , b s := e^ u *' (2t s + e^» x ' y a ) , and Y s :=^^ + y s , 



the following convergences hold, as s — > +oo : 

- . i'2 



6 e -V2u Xs _ i 



2 



f (V.-Y) 



ub. 



2 1 - 



This notion of convergence to the boundary B can be reinforced to a certain extent : 
see Remark [TUl concluding this article. We saw above that any lightlike geodesic belonging 
to a light ray B converges to it. On the contrary, a timelike geodesic does not converge to 
any light ray: by Section [2TTI we get indeed B = (-, |, Y) G B , but the cylinder Cb has 
to small a "radius", since by Equation (8) we have k 2 = [1 — £ 2 — a~ 2 ]/2 < [1 — £ 2 ]/2 . 

Proposition 3 The isometry group of Godel's universe (recall Section^) does operate on 
the boundary B, so that the above definition^ is consistent. Precisely, it works as follows. 

1 ) The translation (t, x, y,z) i— > (t + to, x, y + yo, z + z ) changes (£, g, Y) into (£, g,Y + y ) . 

2) The hyperbolic dilatation (t,x,y,z) i— > (t, x + Xq, ye - ^ wzo , z) changes (£, g,Y) into 

(e,ge y/iuxo ,Ye- y/iuxo ) . 

3) The rotational symmetry (u, r, 0, z) \— > (u, r, + O) z ) changes (£, g, Y) into 

( . ug cos(w0 o ) + [^-a]sin(^0 o ) \ 

IE, a+ lQ-a\cos{ufo) - u gY sin(a;0o), — ] r~T\ • 7 1 m > 

y uj [a + [g — a\ cos(a; O J — ^ g sm{u <po)\J 

where a := 2ach ( 2r )~ s ' 1 ( 2r )^ i s constant under d> i— > d> + <bn , and on each geodesic. We 

u+2sh i r <f> t yj i z 

have indeed : a — |(1 + uj 2 Y 2 ) + , or equivalently : g — a — |(1 — uj 2 Y 2 ) — 1±L _ 



Proof The two first items are straightforward. On the other hand, the action of the ro- 
tational symmetry (it, r, 0, z) t— > (it, r, + <f) , z) is not so obvious. However, a computation 
shows that we have in coordinates (it, r,<ft,z) : a = u + 2(sh r) 2 , and : 

b = A + \P cos(a; 0) — 2c<j~ 1 r sin(u; 0), Z:=cj6F = \E' sin(a; 0) + 2oj~ l r cos(c<j 0) , 

with 

A := 2a ch (2r) - sh 2 (2r) and * := [2a - ch (2r) 0] sh (2r) . 

Note that A = 2a + 4[a — ch 2 r 0] sh 2 r = 2a + 2 J^L is seen to be constant on each geodesic, 
by looking at the expression of the Lagrangian L in coordinates (it, r, 0, z). Alternatively, 
a computation yields : 2A = b + u 2 b (2Yy — y 2 ) + (4a — be"^ 2 ^ 00 ) e _v/2w:c , whence by using 
Proposition [2]: 



2a — g 



2A 



— Q = LU Q 



V 2 _ 2(l-tf 

J .9-9 



4 


to 


1 








+ - 


4- 


1 2 + tg 2 (^)_ 




Q 



4 v / r^i 2 tg 2 (^ 

2 + tg 2 (^) 



4 fl-£ 2 ) 



2 ^ " 2 ' ■> v J ! ; , whence a = | (1 + cu 2 F 2 ) + — 

2 



Now we have at once : under i— > + O > (O; ^; ^) is changed into 

(a , A + [6 — A] cos(co> o ) — Z sin(co> O ) , Z cos(u; 0o) + [b — A] sin(co> O 
so that (£, a, Y) is changed into (recall from the above that a = A/ a) : 

uj g cos(u; 0o) + [a — — ] sin(c<j O 



^ , ^ + [Q ~ „ ] C0S(W O ) - CJ Y Sm(ljJ O ) , 



w [~ + [g - ~ ] cos(a; O ) - w g sin(a; O )] , 



o 



2.3 Ricci curvature and energy tensor 

Recall that the Christoffel symbols are computed by : 



r 



ij 2 



^ ^ 



and that the Ricci tensor (Rij) is computed by 

dT k - r)T k 
d _ »j _ _ ik , r k r e _ r k r e 

n v ~ Qtk Q^j + 1 M L ij 1 if L jk ■ 
From Equations (!'), (2'), (3), we find all non- vanishing Christoffel coefficients 



^xy — Fty — (o> 



2)e^% rl = V2u, r x yy = (u 



2)e 2V2uJX , Y v iT = -^f2ue-^ x 
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Therefore we get all non-vanishing Ricci coefficients : 

Rtt = 2u 2 ; R ty = R yt = 2uj 2 e^ x ; R yy = 2u 2 e 2 ^ x . 

Hence, the scalar curvature is R = g tj Rij = 2oo 2 . 
Einstein equations 

2 9ij ^ 9ij -^i j 

are satisfied, with cosmological constant A = u 2 representing a positive pressure, and 
energy tensor (Tij) = (Rij) = (uiUj), where u := (\^2 u, 0, \[2uj e^ UJX , 0) represents the 
four- velocity of matter, which rotates with constant velocity uj . The energy is thus : 



i + e V2U,X y 



n2 



= 2cu 2 a(£,0 2 . 



3 Relativistic diffusion 



Recall from [F-L J] that the general expression of the relativistic operator £ is : 

£ = + (¥«' - e« J r*(o) £ + 4 («< - /'«)) ^ ■ 

<j being an arbitrary fixed positive (speed or heat) parameter. 

Equivalently in the present setting, the relativistic diffusion (£ s ,£ s ), in coordinates 
£ = (t,x,y,z), solves the following system of stochastic differential equations: 

dt s = i s ds ; dx s = x s ds ; dy s = y s ds ; dz s = z s ds ; 

dt s = -2y/2ui,x a ds- V2iue V2uJXs x s y s ds+ ^i s ds + adM l s ; 
= -V^ue^^isi/sds- (uj/V2)e 2V2u)Xs y 2 s ds + ^x s ds + adM x s ; 
dy s = 2 y/2ue-^ ux ' i s x s ds + ^£ y s ds + a dM y s ; 
rfi s = z s ds + a dM z s ; 

where the M 4 - valued martingale M s := (M*, M x , MJ , M*) has (rank 3) quadratic covari- 
ation matrix : 

i 2 + l i s x s i s y s -2e-^ 

A.2 



( 



(dMldMj) 
ds 



t s X s 



V 



t a ij s - 2 e" 

t<t Zc> 



-V2o 



y s 2 + 2e 



-2\/2 ui x s 



y s z s 



t s z s \ 

X s Z s 

y s is 



Recall that the unit pseudo-norm relation reads : 



(0) 



l + i 2 s +x 2 s + z 2 



il, I 2/', 
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Thus, the relativistic diffusion (£ s ,£ s ) is 7-dimensional, having phase space: 



t,x,y,z,i,x,y,z) G 



1 + i 2 + x 2 + z 2 



£ : = 

or equivalently : 

£ = {(t,x,y,z,i,x } y } z) G 

Note that the particular phase subspace distinguished in Remark [T]: 
£ := £ n [x = y = 0} = £ n [i 2 = 1 + z 2 ; x = y = o} 



y + 2t] } 



l + x 2 + z 2 + U 2 ^ x y 2 



t + e /2wx y 



is clearly not stable under the relativistic diffusion (£ s ,£ s ), contrary to the geodesic flow, 
and even instantly unstable : starting from any point in £q , its exit time from £q is null. 



3.1 Reduction of the dimension 

The study of geodesies induces to consider the following quantities (which, clS Z s , are 
constant along each geodesic), setting (as in Definition [2]) : 

(10) a s := i. + e^ ux 'y a and b s := e^ ux - {2t s + e^ u °» y s ) . 
Then we have : 

da s = ^£a s ds + a dM a s = '^£a s ds + a (dM l s + e^ ux ' dM y s ) ; 

and 

db s = ^fb s ds + a dM b s = '^fa s ds + a e V2u)Xs (2 dM\ + e /2wXs dM y s ) . 
Moreover we have : 

dx s = (u/y/2) e ~ 2V2uJX ° b 2 ds-V2u e" ^ Xs a s b s ds + x s ds + a dM x , 

and the M 4 -valued martingale M s := (Mf, M^, M x , M*) has (rank 3) quadratic covariation 
matrix : 



/ 


a 2 s -l 


a s b s -2e V2ujx ° 


a. 


X s 


a s z s \ 


a £ 


b s -2e V2uXa 


fo2 2 g2\/2uixs 


K 


X s 


b s z s 




(Xg X s 


b s x s 




+ 1 


X s z s 


V 


(X s Zg 


b s z s 


X 


9 %S 


Z 2 s + lJ 



From this, we deduce the following. 

Corollary 1 The (7-dimensional) relativistic diffusion (£ s ,£ s ) admits the following sub- 
diffusions : (a s ) ; (z s ) ; (z s , z s ) ; (a s , z s ) ; (a s , z s , z s ) ; {x s , x s , a s , b s ) . 
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The unit pseudo-norm relation can be written : 

(00) 1 + x 2 + z 2 + (a s - e~^ ux ° b s ) 2 = \ e ~ 2 ^*° b 2 s , 

or equivalently : 

(00') l + x 2 s + z 2 s + \(2a s -e-^ x °b s ) 2 = a 2 s . 

Hence the phase space £ of the relativistic diffusion (£ s , £ s ) can be written equivalently : 
£ = {(t,x,y,z,a,b,x,z) G M 8 I f + x 2 + z 2 + \ (2 a - e ' V2wx bf = a 2 }. 
And the particular phase subspace £q distinguished in Remark [1] can be written: 
£ = £ n [a 2 = I + z 2 ; 2 a = e ~ V2uJX b ; x = o}. 

Remark 4 We see in particular that a 2 > 1 and that b 2 > 2 e 2 ^ U]X ^ for any proper time 
s > . Therefore, (a s ) and (b s ) almost surely never vanish. Moreover, they must have 

the same sign, since (00') implies 
This implies also 



— -2 

a s 



< V2 and then e ~^^ -L > 2 - V2 . 



< 



V2' 



3.2 Study of the one-dimensional sub-diffusions (a s ) and (i s ) 
These two one-dimensional sub-diffusions are easily handled. 

Lemma 1 There exist two standard real Brownian motions (w s ) and (w' s ), and two almost 
surely converging processes (rj s ) and (rj' s ), such that we have: 



exp 



a 2 s + cr w s + r] s for any proper time s > 



and 



exp 



a 2 s + cr w' s + r]' s for any sufficiently large proper time s . 



Proof The stochastic differential equations satisfied by (a s ) and (z s ) are respectively : 



da s = a s ds + o~ \Ja 2 s — 1 dw s , and dz s = ^ z s ds + o J z 2 + 1 (iu^ , 

for two standard real Brownian motions (w s ) and (w' s ). 

These equations are solved as follows ; we have real Brownian motions (W u ) and (WQ 
such that : 



inf jit 


/ 




Jo 
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and 

z s = G(W'[mf {u\ J {I- \W' v \ 2 )- 2 dv > as] 
with F(W) := -S== and G(W) := -JL — . 

Clearly, as u increases to the hitting time of 1 by W, then / (W 2 — l) -2 dv increases 

Jo 

to infinity W[u] goes to 1, and F(W[u]) goes to infinity showing that a s goes almost 
surely to infinity with s. The same reasoning holds for (i s ) (except that \Wq\ must be 
smaller than 1, while \Wq\ must be larger than 1), so that, in the same way, \z s \ goes almost 
surely to infinity with s . The invariant measure of the diffusion (a s ) is l{| a |>i}V a2 — Ida , 
and the invariant measure of the diffusion (z s ) is \/z 2 + 1 dz . 

Then we have almost surely (since \z s \ — ► oo), for any sufficiently large proper time s : 



dlog \a s \ = (1+ 2^j)o" 2 cis+cry 1 — a s 2 dw s , and rflog \z s \ = (l — ^jj)a 2 ds+ayl + z s 2 dw' s . 
Whence, for real Brownian motions w , w' and for sufficiently large proper times Sq, s : 



log \a s \ = log | (So) + cr 2 s + a w s 



and similarly : 



a 



! du 



o 2 a, 2 , 



w 



a a„ du 



log— - = a (s-s ) + a(w s -wj- ir^ + w j— 

\Z S0 \ J s ^Z u LJsol+,/1+2 



o l + Jl-a- 2 



ra a 2 z~ 2 du 



cr 2 s + o(s) > a 2 s/2, 



a 2 s + o(s) > a 2 s/2. 



This implies the convergence of the integrals in the above formulas, and then the result, o 



Corollary 2 For any sufficiently large proper time s we have: \z s 



-o(s 5 / 9 ) 



We have also the following lower control, which we shall use later. 



Lemma 2 For any A > , we have P 
and P[(3s > 0) \a s \ < Ae a2s/2 / \a \ > A 2 



[3s > 0) \z s \ < Ae a2s/2 / \z \ > A 2 ] < l/VX, 
< If A. 



Proof Fix A > v^3 and \zq\ > A 2 . The stochastic differential equation satisfied by 
(log \ z s \), already written in the proof of Lemma [H is equivalent to : 



dlog 



z R e 



-a 2 s/2 



^ (1 - z~ 2 ) ds + a jl + ir 2 dw' . 



Let us apply the comparison theorem (see for example ([I-W], Theorem 4.1)) : setting 



T% := inf{s | \z s \ = Ae a s/2 } and logr s := log A 2 + \ (1 - A~ 2 ) s + a^l + A- 2 w' s 
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we have : 



inf 

0<s<T5 



-a 2 s/2 



> inf r s , whence 



0<s<T 



A 



a<oo]< P[ logr s hits log A] 



P [ w' s - j^=2 s/2 hits log A] = A'^^ 7 < 1/VA. 
Similarly, for |a | > A 2 and T\ := inf{-u | \a u \ = Ae a2s ^ 2 }, since 

-(T 2 s/2 



a., e 



^ a s 2 ds + (T\/l — a, 2 



we get : 



P[T^ < oo] < P[ w s - s/2 hits log A] = 1/A. o 



3.3 Study of the two-dimensional sub-diffusion (a s , z s ) 

Recall from Section 13.11 that we have : a 2 — 1 > z 2 , 

da s = ^Y~a s ds + a dM a s and dz s = ^fz s ds + o dM z s , 
where the Revalued martingale M s := (M",M*) has quadratic covariation matrix: 



a s z s z 2 s + 1 J ' 

We get easily an asymptotic variable of the sub-diffusion (a s ,z s ). 

Proposition 4 The process (z s /a s ) converges almost surely, toward some random limit 
i such that < \£\ < 1 . 

Proof Using Ito's Formula and the above expressions for da s and dz s , we get : 



dz s 
as 



z s da s z s (da s ) _ (da s , dz s ) 



o~ r 
a 



dM* - — dM° 

as 



(J 2 Zs 



ds . 



with 



dMt 



dM? 



\z s /a s 



ds . 



Hence, we have some real Brownian motion W such that almost surely, for any s > 



a 



a 



^rdu + o-W 



o at 



I z u j a u | 



du 



which almost surely converges toward some random limit I e M, by Lemma [H The unit 
pseudo-norm relation (00') implies a 2 — 1 > z 2 , and then £ 2 < 1. 

Similarly, using Lemma [1] again, for any large enough proper time s we have : 



d 



da s a s dz s a s (dz s ) {da s ,dz s } 



o r 



dM a s - — dM z s 

Zs 



+ 



a 2 a s , 
,„ ds , 

zi 
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with 

\a s /z s \ 2 - 1 



i: 1 



dM* - — dMt 

i. 



ds , 



whence the almost sure convergence of (a s /z s ), which proves that £ ^ almost surely, o 

The following statement ensures that the range of possible limits £ in Proposition HI 
is the whole [— 1, 0[U]0, 1]. This provides a continuum of non-trivial bounded harmonic 
functions for the relativistic operator C . 

Proposition 5 For any real £q such that < \£q\ < 1, and for any e > , we have 
p r„ . ,. 

|oq| is large enough. 



o — e < t — lim — < £ + e > 1 — e , provided z /ao is close enough from £ and 



Proof Fix A > 9 , |a | > ^4 2 and \zq\ > A 2 , such that Zo/ao is close to £$ (precisely, 
we demand | log(^-)| < A~ 2 ), and consider the event : 

A := { a 2 > 1 + A 2 e a2s and z 2 > A 2 e^ s for all s > }. 

By Lemma [2] we have : P(.4) > 1 — 2/\TA . Now, on A we have : 



00 du f°° du 9,9 f°° du 9,9 

-+/ — < 2a- 2 A' 2 and / — < a' 2 A~ 2 

o at - 1 Jo zi Jo z± 



a 2 -I 



displayed in the proof of Propo- 



Hence, we see from the expression giving log 
sition HI that we have on A : 

\\og{£/£ )\ < 2A- 2 + a max{\W a \ |0 < s < o' 2 A' 2 } . 

Finally, as 

P[tr max{|V^ s | | < s < o~~ 2 A~ 2 } > A~ 1/2 ] < 2F[max{W s | < s < A' 2 } > A~ 1/2 ] 



2¥[\W A -2\ > A- 1/2 ] = 4P[Wi > VA} < e 



A/2 



we obtain : P[ | log(£/£ ) \ < 2 A~ 2 + A~ 1/2 ] > 1 - 2/VA - e~ A/2 . o 



Proposition 6 The process X s := argsh y\J a 2 — z 2 — 1 J solves the following stochastic 
differential equation : d\ s = a dw s + a 2 coth (2A S ) ds , for some standard real Brownian 
motion w , and then is a real diffusion which goes almost surely to infinity as s — ► oo . 
Moreover, the process fj s := X s — o~w s — a 2 s converges almost surely (inM.) as s — > oo . 
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Proof Using Ito's Formula and the expressions for da s and dz s , we get 



d\. 



2a s 



sh (2A a 

a 2 
sh (2A S ) 



2f-a s ds + adM a s 



2z„ 



sh (2A B 



^z s ds + adM z s 



+ 



a 2 s z 2 s ch (2A 



+ 



ch 2 A< sh 2 A 



A — 1) - gh(2As 



+ 



ch 2 A, sh 2 A 



ch 3 A S sh 3 A 
2 - 1 

(zi + 1) rfs 



^ 2 j 

a as 



2a 



sh (2A S ) 



3<7 2 ch 2 A s , a 2 i 2 ch(2A s ) 2 , 



sh (2X S 



ch 3 A, sh 3 A, 



a' 



sh(2A s ) 



a s + z* — ch 2 A S ^ a 2 i 2 + sh 2 A S 



ch 2 A, 



sh 2 A, 



ds 



2a 



sh (2A S 



-((a 2 - i 2 )(a 2 - i 2 - l)) 1/2 dw s + 



sh (2A S 



3ch 2 A, + 2 



a 2 i 2 af + if 



ch 2 A, ch 2 A< 



= cr £#tt) s + a 2 coth (2A S ) ds . 
Since coth (2A S ) > 1 , the comparison theorem ensures that we have almost surely : 

A s > A + o w s + cr 2 s — 



-oo . 



Moreover, we have almost surely for large enough s and for s > Sq: A s > a 2 s/2. 

rs 

Hence, we deduce that fj s = fj sn + 2<r 2 / — rr converges almost surely o 

Js e 4Au — 1 



Remark 5 The equation satisfied by (A s ) can be precisely solved as follows, provided 
A > , using some real Brownian motion (3 started from /3 = \ log[coth A ] : 

A s = \ log coth (j3 inf |w J sh.~ 2 (2/3 v ) dv = a 2 s| 

This implies that we have almost surely : \ s > for any s > : the state subspace £ 
of Remark [1] is polar for the relativistic diffusion. Recall from Section [3] (before Section 
13. ip that it is also instantly unstable. Hence we can always restrict the state space S of 
the relativistic diffusion to S B\ Sq . 



By the unit pseudo-norm relation (00') and by Remark [U we have almost surely for 
any s > 0: a 2 — 1 > z 2 . Therefore there exist two independent standard real Brownian 
motions (w s ) and (w s ) such that : 



da x = '^y~ a s ds+a\/a 2 — 1 dw s and dz. 



z s ds+a 



d 2 - 1 



dw s + a 



\ 



z 2 

1 5 - ? — - 

a 2 - 1 
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We shall need to know that in fact \l\ < 1 almost surely. Let us set 



11 



\ 



1 



1 + z 2 





2 + l 

~ 2 


r b - 

2 — e -V2ux s _£. 


2 


sh A s 






L a s - 







Note that the phase subspace Sq (of Remarks [T] and E]) is precisely : £q — £ H {A = 0}. 

Proposition 7 T7ie random limit £ = lim s _ 4 , 0O (i s /a s ) of Proposition^ satisfies almost 
surely: < |£| < 1 . 

Proof By Proposition H] and Notation (11), we have almost surely : log(l — £ 2 ) = 2 log-A,^ . 
On the other hand, using Ito's Formula and Proposition [6j we get : 



d(\ogA s ) = d(log[sh A s ]) — d(log|a s |) = cothX s d\ s 



d(X s ) 
2sh 2 A< 



-a 2 {l + \a- 2 )ds-—dMS 

Ojc 



a 



sh 2 A s a. 



'" )dM? --^—dMt 



sh 2 A< 



+ 



2 coth A s coth (2A £ 



sh 2 A s 



-2- 



ds 



a 



a s sh 2 A 

whence, for some real Brownian motion B : 



z 2 + l) dM? - a s z s dM* 



——t ds , 

2 a 2 ' 



lo K (l- 



2 log A + 2a B 



z 2 + l 



o at (a 2 — zj — 1) 



ds 



a 



00 ds 

9 ' 

o a. 



which converges (in K) almost surely, by Lemma [T], Proposition H] and Proposition 0, 
showing that indeed £ 2 < 1 almost surely, o 

We have furthermore the following. 
Proposition 8 The law of the random limit £ = lim (z s /a s ) has no atom. 

Proof Fix any £ £ ] — 1, 1 [ , and set 6 S := z s — £oa s . The stochastic differential equation 
satisfied by (S s ) is easily seen to be : 



d<L 



5 S ds + a J5 2 +l-£ 2 d(3 s 



for some standard real Brownian motion (/3 S ). This diffusion equation can be solved as 
follows : we have a real Brownian motion (W u ) (started from Wo £ sucn that : 



inf jii 


/ 




Jo 



r« (1 - If) dv 
-(l-t 2 o) 2 W v 2 



> a s 



}]) 
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with F(W) 



then 



{1 _f2 )3 /2 w 



As u increases to the hitting time of ±(1— £q) 1 by W, 



increases to infinity, W u goes to ±(1 



and F(W U ) goes 



lo \-{\-£\) 2 W 2 

to ±oo , showing that \b~ s \ goes almost surely to infinity with s . (The invariant measure 



of the diffusion (5 S ) is y/6 2 + 1 - £\ dS .) 

Then we have almost surely, for any sufficiently large proper time s : 

d\og\5 s \ 



1 - ^)o 2 ds±o\ll 



-gp- d(3 s ■ 

Whence, for real Brownian motions w, w and for sufficiently large proper times Sq, s : 



log 



\5„ 



a 2 (s-s )+a (w s -w SQ )- 



^ 2 (i-^) 



so 0., 



du 

0,. 



a 2 s+o(s) > a 2 s/2. 



This implies the convergence of the integrals in the above formula, and then the existence 
of a standard real Brownian motion (w e s °) and of an almost surely converging process 
such that almost surely, for any sufficiently large proper time s we have : 



to a s 



exp 



a 2 s + aw e s ° +r] e s ° 



exp 



a 2 s + o(s 5 / 9 ) 



For the same £0 and (5 S ) as above, we get as in the proof of Proposition 0]: 



log 



\S s /a s 
\S S0 /a s 



-a 



z u du a z 



o~u 



1 - (z u /a u ) 2 du 
51 



+ aW 



1 - (z u /a u ) 2 du 
51 



almost surely for any sufficiently large so, s . Using the above, this shows the almost sure 



convergence of log 



hence by Proposition 01 that indeed ¥[£ 



0. o 



The following statement implies that the asymptotic a-algebra of (a s ,z s ) is generated 
by the only variable £ . 

Proposition 9 The process (z s — £a s ) converges in law, toward some smooth law, but not 
in probability. Furthermore, we have (z s — £a s ) = o(\a s \ 2 ^ 9 ) almost surely. 



Proof Since 



dM? 



[1 - a s )ds 



dM z s - — dM a s 
a 



1 



ds , and 



dM« z s ,.,„\ z. , 

dM* - — dMf) = ^ds, 



a q a q ' at 



there exist two independent standard real Brownian motions w, w' such that : 
dM* 



1 - a~ 2 dw s , and dM* - — dM* 

a. 



a1\ 1 — a: 



dw s + 



at — 1 — z 



>2 



\ <*2 



dw' . 
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Then, since 



dloga s = a 2 {l + ^n)ds + —dM^, 

a., 

we have for any < s < u (recall Lemma [T]) : 

<' s + u dv 



cxp 



— a 



s+u 



-o- 2 «-o-(io s+u -«) s )+o(l) 



Hence, the expression for c?[i: s /a s ] used for Proposition HI implies that almost surely: 



Zc Z (La 



a 2 a s / du — a a s 
Js at 



dM> - ^ dM a u 



a 2 a s f 0(a M 2 ) du — a a s ( Wu — a a s f 



1 - 



dw[, 



\ al - 1 a u 



Now, there exist Brownian Motions W, W (independent of s) such that : 



w 



alJl - a-~ 

(I- 



z 2 du 



00 z\ du 



«5(1 - O 



4/9 



o a 



|- lfl /®) X 


" /"OO 


a s 


4 


4 / 9 

= o(|o,r i6 /») x 


[/ 




Jo 








Jo 



4/9 



4/9 



o(|a s |- 4/3 ) x f°° e-^-^^-^^'-^du 



4/9 



(|a s |- 4 / 3 )x(^exp(-[7^+( S + W )(l + f^) + S (l-^)]^V2)^) 4/9 = o(|a. 



= o[\a s 
and similarly : 



-4/3 



zl dw' u 



\ a\-\ a u 



W' 



z 2 , i du 



al - 1 



a u du 



4/9 



= o(\a s 

(l„ |- 2 /3 



-8/9 



X 




a s 


2 

du 


4/9 

= o(|a,r 8 ' s ) X 


[/ 




Jo 








Jo 



( g-2o- 2 w-2o-(«) s+u -io s ) ^ u 

Jo 



4/9 



o a 



x [ I exp 



7 u+ (s+u) 14 



a 2 /4Adu 



4/9 



o a. 



1-7/9 



This shows also that J @{a u 2 ) du = o(\a s \ 7//4 ). 

So far, we have shown that almost surely, as s — > oo : 



z.-ta, = o a 



a a. 



<k| 2/9 )- 
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Considering then a standard real Brownian motion to", independent from (a s ,z s ), we 
have : 

, r°° n r I ~ " 

dw" 



s a u 



-2cr 2 M-2o-(ui s+u -ui s )+o s (l) 



-vT 



1 2 



< 



-2<T 2 « 



duX J J e -2a 2 u-4a(w a+u -w s )+o s (l) ^ 



which goes to in probability as s — > oo . Hence, (z s — £a s ) behaves in probability as: 

/oo „ 

and then converges in law, toward the law of : 



roc / roc 

VT^¥ J e~ u - w ^dwl = Ml ~P)J e~ 2u ~ 2Wu du x N , 

N denoting a A/"(0, 1) Gaussian variable, independent from (i,w). 

On the other hand, the above proves also that, as t > s — > oo , (z t — £a t ) — (z s — £a s 
behaves in probability as : 



VT 



a 2 {u-s)-o(w u -w s ) ^ w ll _ l e -cr 2 (u-t)-a(w u -w t ) ^yjtt 



which diverges in probability : for small enough e > and large enough (t — s), we have 



P 



-a 2 {u-s)-a{w u -w s ) ^Jl 



< e 



P 



> 1 — e , and by independence : 

> 2e 



OO 



e ~cr 2 (u-s)~cr(wu-w s ) ^11 _ I ' ^a 2 (u-t)-a(w u -wt) 

U Jt U 



> 1 - e . o 



Propositions El E] and Notation (11) imply at once the following. 



Corollary 3 We have almost surely : A s = \J 1 — £ 2 + o(|( 



I -7/91 



3.4 Study of the four-dimensional sub-diffusion (x s , x s , a s , b s ) 
Recall from Section 13.11 that we have db s = b s ds + a dM h s , and 



dx s = —( e -^" Xs — - 2)e-^ ux 'a a b s ds + ^x s ds + a dM x 
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the R 3 -valued martingale M s := (M£,M%,M£) having quadratic covariation matrix: 



a s o s — ze v s a s x s 
a s b s -2 x ° b 2 s - 2 e 2V2u Xa b s x s 

b s x s xl + l 



QjQ Xc 



The process (b s ) alone is easily handled, analogously to Lemma El 

Lemma 3 There exists a real standard real Brownian motion (w"), and an almost surely 
converging process [f]"), such that we have: 



cxp 



o 2 s + a w" s + rj" for any proper time s . 



Proof We already noticed in Remark H] that the unit pseudo-norm relation (00') forbids 
any vanishing of b s . Thus there exists a real standard real Brownian motion (w") such 
that for any proper time s we have : 



d log 1 6,| = a 2 ds + a 2 e 2V2uJXs b~ 2 ds + a y/l - 2 e 2 ^^ b~ 2 dw" s 

so that there exists another real Brownian motion w" such that : 

du - f <-s[ ^V2, 



lo 



g \b s \=\og\b \+a 2 s + a 2 f S e 2 ^ x ^ + aw': + V2aw"l f 

Jo b\ \Jo 

Now using Remark H] and Lemma Q] we get : 



b 2 , + \ b 2 -2e 2 ^*^ 



du] 



bl 



poo 




lo 


^ 6 J 



a.. 2 du < oo 



which shows that rj" := log |6 S | — o~ 2 s — o w" almost surely converges as s — > oo . o 
Then we get easily a new asymptotic variable of the relativistic diffusion. 

Lemma 4 The process log(b s /a s ) converges almost surely as s — > oo . 

Proof Recalling from Remark [4] that b s /a s > 0, we have for any proper time s > : 



dlog 



a 2 e 2V2wx ° b~ 2 ds-\a 2 a" 2 ds + a (b; l dM b s - a~ l dM i 



or equivalently, for some real Brownian motion W : 



\ h JL 


-log 


r6o- 


= * 2 / 






-a - 


■/ 



6 2 



o 2 a,, v ./o 



Qi U b b 



b 2 



du^j. 
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Now, as we already noticed in the proof of Lemma [3] for / e 2 ^ 2wx ^ by Remark [4] and 

Jo K 

Lemma [1] we get : 



(Xti. b lt 



du 



poo 




Jo 


t bj 



a,, 2 du < oo 



and then the almost sure convergence of log(o s /a s ) as s — > oo . o 
By Remark H] again, we deduce at once the following. 

Corollary 4 The process (x s ) is almost surely bounded. Setting g := lim — , we have 
precisely : 



< (l - -^=) q < limmf e^ ux - < lim sup e^ wx ' < (l + -^=) g < 



oo . 



The following statement, analogous to Proposition \5[ ensures that the range of possi- 
ble limits q in Corollary H] (and Lemma H]), is the whole ]0, oo[. This provides another 
continuum of non-trivial bounded harmonic functions for the relativistic operator £, . 



Proposition 10 For any real Qq such that < g < oo , and for any e > , we have 

r b s i 

F po — £ < Q = lim — < Qo + s > 1 — s , provided b /a is close enough from g and 
1 s ^°° a s J 

| do | is large enough. 

Proof Fix A > 1 , | ao | > A 2 , and yo such that Yq is close to go ; and use the expression 
displayed for log(b s /a s ) in the proof of Lemma HI Remark HI and Lemma [2j to get on an 
event of probability > 1 — 1/A: 



log q - log 



a 



< ° 2 d + in? 



du 



o at 



+ a max \ \ W u 



0<«<4(1 + ^) 



du 



o at 



< 3 A' 2 + a max{|W u | |0 < u < 7 a~ 2 A~ 2 } , 
b 

so that P( | log q - log g I < log f— — 1 + 3 A~ 2 + A~ 1/2 ) > 1-1/ A 



-do Qo 

Recall that we set (in Section 13.31) : 



e- A / 14 . o 



(11) A s 
Let us set also : 



x s 



+ 



2 - e~^ u 



\ 



1 + z 2 



sh Ac 



(12) x s = a s A s cos 7 S ; 2 — e 



, b s 



V2A S si 



sin 7 S 
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Remark 6 By Corollary [3X we have almost surely : 



e -V5 w *. _£ = 2 _ J 2 (l-p) sin 7s + o(|a s |" 7/9 ) and x, = v 7 ! 3 ^ Gg COS7s + o(\a s \ 2 / 9 ) . 



Proposition 11 There exists a standard real Brownian motion W such that (provided 
Aq > 0) we have almost surely, for any s > : 



ls = 7q + U [ e ~V2u;x u budu + aW 

Jo 

Proof Let us differentiate the relation 



du 



cotg7 s = V2 



X u 



2-e 



0.s 



which makes sense almost surely at least for a dense subset of positive s . We get : 



(1 + cotg 2 7 s )(-<i7 s + cotg7 s (dj s )) = \Pl 
+ V2 



^d 

a s 


^-\f2ux s bs_ 
a s _ 


+ (d 


is 
a s _ 


,d 


e ~V2u)X s bs_ 
a s _ 


) 




2 — g-V^wXj ks_ 
a s 


2 



2 _ e ~V2i0x s ^f 1 ' 



a s 



2 — g"v/2<*>iE8 bs_ 



3 > 



with 



X s 



doc § Xj § dcL § \ d/Q/ g . doc § ) »3j ^ t dcti g / 
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g-^w*. — — 2 



a, 



3> c 



dMf - — dM s a 
a. 



and 
d 



~ e ~V2wx s ^± 









2a 2 



, a 

ds-\ e 



-y/2<jj x a 



dM b s - — dM a s 
a. s 



Hence we get 
V2A 2 S (rf 7s -cotg 7s (d Ts )) 



6. 



- \/2 u) x s " s 2 



— -2 



\/2 a* 



'/- 



+ 



— -2 



(7 



dMf - — dM a s 



2(t2 ^ _ e -V2cxsfMl _ ^ we -Vi«x.^ 



ax 



dM° s - — dM a s 
a* 



2e 



\/2(, 



Xe 6e 



a, at 



2-e 



-1 O X s 2^/2 uxa 



h h 2 



a, at 



'/> 
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a 



dM x -—e-^ x ...... 

a„ 



X u 



dM° + 2 — dM" 



+ V2oje~^ Xs b s A 2 S ds 



+ 



-\/2ujx 3 °« 



— -2 



o 2 x. 



■ds + 



-V2ujx a °s 



— -2 



-i a 2 i s 



I, 1,2 
^ e -V2uix a _ e -2 v / 2a;z s 1£ _ 2 



This implies 



o~ 



2At(d ls ) = - 



-V2ux 3 bs 



— -2 



dMT - — e -^ Xa dM\ + 2 — dM«) = 2a 2 a7 2 A 



2 -2 ,i2 



"8 S J 



whence for some standard real Brownian motion W : 



d'js 



a 



A. a, 



■dW. + ue-^ ux ' b s ds + 



a. 



a 2 x s 



V2w x 3 ^ _ 2" 


-1 cr 2 x s 


2 - 




y/2A*<fi. 





h 

^ e -V2ujx a ^s_ ^ _v v , , 



-2^2 



V2Alal 
b 2 



ds 



ds 



a 



A* a. 



■dW.+ue-^" 3 * b s ds. 



Recall from Remark [5] that (provided A > ) A s a s almost surely never vanishes. This 
shows that the above formula almost surely holds for any s > , proving the result. 

By Proposition [6] and by the definition (11) of A s , we have almost surely : 
A 2 u a 2 u = sh 2 X u = sh 2 [a 2 ?i + aw u + fj u ] , 

/CO 
— — - . Hence we deduce the following. 

Corollary 5 There exists a converging process f] such that we have almost surely, for 
any s > : 

Is = oo I" e'^^budu + fis. 
Jo 

Using Corollary Corollary HJ Lemma [3] and Remark [HI we deduce also easily the 
following. 



Corollary 6 We have almost surely, for large s : 7 S = sign(&o) e 



CT 2 s+o(s 5/9) 



-V2ux 3 = 2 V^-t 2 ) d / c a 2 s+o(sh 



Q sign(b )£> 



(X - 



In the same vein as Proposition [9], we have the following. 
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Proposition 12 We have (b s — ga s ) = o(\a s \ 2 ^ 9 ), almost surely as s — > oo 



Proof We have : d 



2a z 



a 2 bs A i ° 
ds — a — as H 



dM b s - — dM a 
a s 



dM b s - ^ dM s a \ = (4e^ wa *^ - 2e 2V3 
a s / v a s 



-|)ds 



dM? 



[1 — a s 2 )ds . 



and 



dM a s , dM s b - — dM s a ) = (— - 2 e^^ds . 
a s a s 



Hence there exist two independent standard real Brownian motions ty,ii/ such that 



dM„ 6 - — dM? = , dw s + 



at a/1 — a: 



dw' 



and as in the proof of Proposition 0: 



dM s ° r 

= y 1 — a s z dw s , and 

a<j 



-a 2 u-cr(w s+u -w a )+0(a s X ) 



Therefore, we get almost surely : 



q = 2<t 2 / — du — er 2 / -| du + a 

a s a,, Js a,, 



dm,, 



al a/1 - a; 



2 (^ 2 + ^ 2 ) , / 

2 1 

af, - 1 



Now, as in the proof of Proposition [9j we have 



,\/2(jjx u 



du 



f°° b t \ r c 

+ / du = o(\a s \~ 7 ' 4 ) , and / 
Js at v ' Js 



b u -2a u e ^ x ™ 



a%\/l — a., 2 



dw u = o(\a s \ 4/3 ), 



and 



2(zl + xl 



ft dw' u = o(\a s \- 7/9 ) ■ 



So far, we have shown that almost surely, as s — > oo : 



b s - ga s = o(\a s \ 1/3 ) - era 



2 (i^ + £2 J , . v 
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Remark 7 Almost the same proof as above shows that we have in fact (b s — ga s ) = 
o(\a s \ e ) almost surely, and similarly, (z s — £a s ) = o(\a s \ e ) almost surely, for any e > 0. 
But Propositions l9l [T2l (as Corollaries [2, O [6]) as stated are sufficient for our purpose. 

It is possible to complete Proposition [T5] in a way somewhat similar to Proposition M'. 
starting from the last formula above (in the proof of Proposition [T2"l) . and proceeding as in 
the proof of Proposition [9], we can see that the process (b s — ga s ) asymptotically behaves 



as 



a a. 



lu \ 



-/2 I rtf 2 



dw" 



for a standard real Brownian motion w", independent from the filtration of (a s , b s ,x s , z t 
and then, using the pseudo-norm relation (00') together with Proposition [121 as 

r 2„2 



era. 



gV2 uxu 2 g2v2wa;u q2 



dw" 



W" 



Age 



y/2u 



— 2e 



era: 



du 



Note that to complete the Brownian representation of the martingale (M s ) in the proof 
of Proposition [T21 we can find a Brownian motion w'" in the filtration of (a s ,b s ,x s , z s ), 
independent from w, w', such that 



dM? 



Cl s X s 

al-l 



(b s -2a s e V2wx - 



X c 



2_x J2e 2 ^»*s(z 2 e +x 2 e ) 



dw' + 



zi + ±1 



: dw"' . 



Now 



Jo a z s+u Jo 



'O . — /-i 2 

= 2V2w(a; s +u-a; s ) |fa ^ 



2 j i 2 ' 



converges in law to 



() , g\/2w x u -2(t 2 u-2<tw u ^ u f ^2\f2uj x u -2a 2 u-2aw u 

o ' io 



dw , / e 

o 



-2a u—2aw n 



du, V 



where V 1 := - — V 2 ^ 1 e ) s { n jj ; [/ denoting an independent variable, uniform on the 
circle. We get thus the convergence in law, but not in probability, of (b s — ga s ), to the 
smooth : 



W" 



a 2 r kg Ve^ Xu - 2 y 2 e 2V2ulXu - g z 
Jo 1 



-2a u— 2aw u 



du 



3.5 Irreducibility 

Proposition 13 (i) The relativistic diffusion is irreducible : from any starting point, it 
hits any non-empty open subset of the phase-space £ fi\ £ with a strictly positive probabil- 
ity. 

(ii) For any starting point (in £), the law of the asymptotic variable (£, g) charges any 
non-empty open subset of the range (] — 1, 0[ U ]0, 1[) x]0, oo[ . 
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Proof (i) We know from Proposition [T] (in Section 12. ip that there are piece- wise geodesic 
timelike continuous paths, and then trajectories in the support of the relativistic diffusion 
(£.,£.), moving at will the coordinates (t,x,y,z). 

Owing to the the quadratic covariation (rank 3) matrix of the R 3 -valued martingale 
(Mf,Mg,Mg) (recall Section I3TT1) and to the unit pseudo-norm relation (00), we can find 
three independent standard real Brownian motions (w 1 , w 2 , u; 3 ) such that: 



dz s = 3 -|- z s ds + a J z 2 + 1 dw\ 



da s = d s ds + a Qs Zs dwl + a 



*7i • \ ^ 2 +i 



"~ : 1 dw 2 . 



dbs = ^ bsds + a dw i + (j « A^2e^E ±l) dw 2 + a V2e^»*«x. 3 

v/(i 2 + l)(a 2 -i 2 -l) " ^ - i 2 - 1 

Let us use the support theorem of Stroock and Varadhan (see for example ([I-W], 
Theorem VI. 8.1)). We see thus from the above stochastic differential system, that the 
following trajectories belong to the support of (£.,£.) = (t.,x., y.,z.,z.,a., b.,x.) : 

- trajectories moving at will the coordinate z , without changing the coordinates (t, x, y, z) ; 

- trajectories moving at will the coordinate a , without changing the coordinates 

(t,x,y,z,z) ; 

- trajectories moving at will the coordinate b , provided z ^ , without changing the 
coordinates (t, x, y, z, x, a). 

So far, it has become clear that it is possible, within the support of the relativis- 
tic diffusion, to move any point of the phase space £ having given first coordinates 
(t, x,y, z, z, a) G M 6 , onto some point of the phase space £ having prescribed first co- 
ordinates (£', x',y',z',z',a!) G 1R 6 . 

It remains only to consider the last two coordinates (b, A). They are of course con- 
strained by the unit pseudo-norm relation (00'), which tells precisely that they run some 
ellipse of this plane of coordinates, which is centred on the axis {A = 0}. The last type of 
trajectory mentioned above allows now to move (b, A) arbitrarily on the upper half and on 
the lower half this ellipse, without changing the other coordinates, within the support of 
the relativistic diffusion. 

Finally, we must make clear that we can cross the axis {A = 0}, within the support 
of the relativistic diffusion, at the cost of an arbitrarily small move of all coordinates. 
Now, owing to the the quadratic covariation (rank 3) matrix of the R 3 -valued martingale 
(M",Mg, Mf) (recall Section I3TT]) and to the unit pseudo-norm relation (00), we can find 
three independent standard real Brownian motions [w l ,w 2 ,w z ) such that: 



dx s = (u/y/2 ) e- 2V2uj Xs b 2 s ds-V2u e _v ^ wx * a s b s ds + ^ x s ds + a Jx 2 s + 1 dw] ; 
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da* 



db s = ^Y~b s ds + a 



o 2 7 ^ S i 1 

^p- a s ds + a dw s + a 



he Xa 



±i + 1 



: dw s + a 



a s b* — 2e 



±1 + 1) 



^ + 1)^-^-1) 



dw„ + cr 



a? — xl — 1 



This shows that, arrived in an arbitrarily thin ^-neighbourhood of the axis {x = 0}, we 
can cross this axis (within the support of the relativistic diffusion, acting on the Brownian 
component w 1 ) without changing the coordinates (t, x, y, z), and perturbing the coordinates 
(a, b, x, z) only by some move of order 5 . This ends the proof of irreducibility. 

(ii) This is a direct consequence of (i) above and of Propositions 151 and [TOl: by (i), 
it is indeed enough to start the relativistic diffusion so that z /a Q be close to a given 
£ o £ (J — 1,0[U]0, 1[), b /a be close to a given g > , and |a | be large enough, o 

3.6 Convergence to a lightlike geodesic 



From Section [2721 (Proposition [2] and Remark [3]), for a lightlike geodesic we have three 
geometrically meaningful real parameters £, g, Y such that : 



— = £e[0,l]; - = f?e]0,oo[; ^-x T + y T 
a a uj b 



Y e 



and 



6 -\/2ujx t 



+ 



uj g 



(Vr - Y) 



where r denotes an (irrelevant) affine parameter. 

Accordingly, let us exhibit now a third asymptotic random variable Y for the relativistic 
diffusion. 

r~ X 

Proposition 14 The process Y s := Us + ^ t~~ converges almost surely, as s — > oo ; 
toward some real random variable Y . And we have Y — Y s = o(|a s | _4//9 ) almost surely. 



Proof Recall from Formulas (10) that we have y s = e v ^ UJXs (2a s — e ^ UXa b s ). We have 
then : 



d 



dx s x s db s x s (db s ) (db S} dx s ) 
be b 2 „ b\ b 2 „ 



j, e -2V2u JX3bsds _ y /2 ue -V2u JXs ^ ds _ 2a 2 e 2V2u Xa ±s ^ + O _ a tl dM 

Vi M o„ bi 



S 1 



whence 



V2 



dY s 



_ 2a 2 e 2V2. Xs ±± ds+ ° dM *_ a ±± dM b 

65 6- bl 



29 



and for some Brownian motion W 

■2 



V2 



V2 



2-/2wx„. X u 







bt 



/* 






2- 


du 


Jo 




b u 




u . 



By Corollary HI Remark and Proposition [T2J the two above integrals, and then Y s 
converge almost surely, and moreover, we have (for some Brownian motion W) : 



Y - Y = Y - Y 



2^/2c 



30 e 2V2^u^ du + W 



6 3 



2 a 2 


\ — ^\/2u)Xu 




2- 


du 


_^ Js 




b u 







0(i; 



a~ 2 ds + W 



0(1)/ a~ 2 du 



a„ 2 du 



4/9 



(«; 8/9 ) 







a s 


2 


) 






du 




Jo 


O's+u 





4/9 



(a; 8/9 ) 



2a 2 u-2a(w s+u -w 3 ) j 



a 



4/9 



o («; 4/9 ) 
(«; 4/9 ) • o 



'du 



4/9 



Proposition [TH together with Remark [6] and Corollary HI implies easily the following. 
Corollary 7 We have almost surely, as s — > oo : 



6 --J2wx s _ i 

2 



+ 



(v. - 



The following statement, analogous to Propositions 151 and [TTJ| ensures that the range of 
possible limits Y in Proposition [bil is the whole R. This provides again another continuum 
of non-trivial bounded harmonic functions for the relativistic operator C, . 

Proposition 15 For any real y and any e > 0, we have F[y — e < Y < y + e] > 1 — e , 

provided Y is close enough from y and \qq\ is large enough. 



Proof Recall from Lemma [2] that the event A := ||a s | > \J\clq \ e a2s ^ 2 for any s > oj has 
(for | a-o | > 3) probability larger than 1 — |ao| -1 ^ 2 - The proof of Proposition [bil shows that 



Y-Y \ = 0(1) f -£+wax{\W.\\0<8< 0(1) f -£} = 0(^) + W{fc) on A, 

> that F(\Y - Y | < 2 |a |~ 1/3 ) > 1 - 2 |a |~ 1/2 , for large enough \a \ . o 
Proposition [TBI improves Proposition (|T3| (n')). We deduce indeed at once the following. 



Corollary 8 For any starting point (in £), the law of the asymptotic variable (£, g,Y) 
charges any non-empty open subset of the range (] — 1, 0[u]0, 1[) x]0, oo[xR . 

More precisely, if the starting point of the relativistic diffusion satisfies: io/ a o close 
enough to £o e] — 1, 1[ ; b /a close enough to g > , Y close enough to y 6l , and \a Q \ 
large enough, then with arbitrary large probability, (£, g,Y) is arbitrary close to (£o,Qo,y)- 
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Remark 8 A rapid look at Remark [3] could let think that there could be a fourth asymp- 
totic random variable for the relativistic diffusion, namely a possible almost sure limit for 
X a := z s + £t s — {f-/^)ls • But as a matter of fact, there is no such limit, in accordance 
with the last sentence of Remark [21 on the geometric irrelevance of additional parameter 



[Z ,T ). Indeed, since by Equations (10) we have i s 



V2t 



a, , we deduce from 



Corollary that X s 



' a u )du converges almost surely, and then, by (the proof 



of) Proposition [9j that so does also 

P S P oo 

X s + aVT^¥ / / e-^M-^-"-) dw" v du . 

JO Ju 

Now, as (setting wr> := w n+ . — w n , for any fixed n e N) 



n+l /-oo 



-° 2 (v-u)-v(w v -w u ) d ,f du = w n 



min{l,i)} 



du 



e- 2 ° 2v - 2 ™« n) dv 



has a constant law, we see that (X B ) indeed cannot converge in probability. 



The theorem of the introduction (section [T]) is now established. Indeed, gathering 
successively Remark [5] and Proposition [131 Propositions H] and d Corollary HI Proposition 
H4l and Corollary [71 and Corollary [HI we get the following main result (for which a > is 
necessary, due to the observation made after Definition [21 in Section 12. 2j) . 



Theorem 1 (i) The relativistic diffusion is irreducible, on its phase space S fi\ So ■ 

(ii) Almost surely, the relativistic diffusion path possesses a 3- dimensional asymptotic 
random variable (£, g, Y), and converges to the light ray B = (£, g,Y) 6 B , in the sense of 
Definition [H Indeed, we have almost surely, as proper time s — ► oo : 



z s /a s 



G]-1,0[U]0,1[; b s /a s 



2 



- 1 



+ 



g e]0,oo[ ; Y s - 
- J(l-0- 



Y e 



(Hi) The asymptotic random variable (£, g, Y) can be arbitrary close to any given 
(£o, Qo,y) ^ ] — 1) l[ x ]0, oo[xl, with positive probability. Hence, the whole boundary (space 
of light rays) B is the support of light rays the relativistic diffusion can converge to. 



Remark 9 From the proofs of Lemma HI Proposition [HI and Proposition HI we have the 
following representation of the asymptotic variable B = (£, g, Y) : 



„ du 



POO _ 

, 2a 2 / — -a , 

a Q Jo az Jo a 6 



o i°° b u du r 00 ! 

Jo 



dM h u - ^ dM« 



Y = Y n 



o of, w Jo 



dMl - % dMl 
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dM" - ^ dM" . 

By Proposition [BJ the law of the asymptotic variable B has no atom, and by Theorem 
(JTJ(m)), it is really three-dimensional. None of £, g,Y is a function of the two others. 

Recall that the random excitement of the relativistic diffusion is a standard three- 
dimensional Brownian motion. Therefore, Theorem HJ reinforced by Proposition [9] and 
Remarks H, [HJ EH incites to believe in the following. 

Conjecture The tail cr-field and the invariant a-field of the relativistic diffusion in Godel's 
universe are the cx-field generated by the asymptotic three-dimensional random variable 
B = (£, Q,Y) of Theorem [T] (exhibited by Proposition HJ Corollary SI and Proposition UM . 



£ = a 2 -^du + a a u l 



Remark 10 The convergence of Theorem (Pfl(ii)), of the generic diffusion path (£ s ) to 
some light ray £ = (t,x,y,z), in the sense of Definitional occurs in fact in some stronger 
sense. Indeed, by Remark [6], Propositions [12] and HJ and Lemma [IJ we have on one hand : 



-V2^ s = 2 _ V^Efl sin7s + o( e - 2 ° 2s / 3 ) and y s = Y 



LU Q 



cos 7 S + o(e 



-2a 2 s/3 



while by Remark [31 using the increasing diffeomorphism (p = (r i— > (p T ), we have on the 
other hand : 



2 y/^P) ■ 



sin 7 S and y 



Y 



Q Q " *<P 1 (7s/") " U1Q 



Hence, we have in the (x, y)-plane a strong convergence, of the projection of the generic 
relativistic diffusion path to the projection of a lightlike geodesic, in the Skorohod topology : 



+ 



o( e 



-2a A s 



/»). 



Otherwise, by Proposition [91 Corollaries [5] and [61 and Remark [3X we have : 

z s +£t s = z s +£t +£ (\e~^ x -b u -a u )du = A 7s + 0(1)+ (\z u -£a u )du = i ls + o{\ ls \ 1/z ) 



z<p-H7.M + fty-Hy./!*) + °( W 1/3 ) = + ^ V l (7./<*)) [ X + °( e CT2s/2 ) 



oo . 



Hence, again in the Skorohod topology, and in the (z, t)-plane, the projection of the limiting 
light ray stands for an asymptotic direction for the projection of the generic relativistic 
diffusion path, but there is no exactly asymptotic lightlike geodesic (only a parabolic branch 
occurs). 
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